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ABSTRACT 

Multifractality of charged particle pseudorapidity distributions is analyzed in central collisions 
of carbon and copper nuclei at 4.5 GeV/c per nucleon. Within the method of normalized factorial 
moments, modified to remove the bias of infinite statistics in the normalization, intermittency-like 
behavior (scaling) is observed up to eight-particle fluctuations. The study indicates a possible 
non-thermal phase transition inside cascading and two different regimes presented in multiparticle 
production. Dynamics of fractal structure formation is further studied using the chaoticity approach 
recently proposed. The distributions of the horizontal factorial moments are considered and a scaling 
behavior, referred to as erraticity, of the normalized moments of the distributions is obtained. The 
corresponding entropy indices are calculated indicating chaotic nature of multiparticle production 
with a specific self-similar structure. 
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Multifractality of charged particle pseudorapidity distributions is analyzed in cen- 
tral collisions of carbon and copper nuclei at 4.5 GeV/c per nucleon. Within the 
method of normalized factorial moments, modified to remove the bias of infinite 
statistics in the normalization, intermittcncy-like behavior (scaling) is observed up 
to eight-particle fluctuations. The study indicates a possible non-thermal phase 
transition inside cascading and two different regimes presented in multiparticle 
production. Dynamics of fractal structure formation is further studied using the 
chaoticity approach recently proposed. The distributions of the horizontal facto- 
rial moments are considered and a scaling behavior, referred to as erraticity, of the 
normalized moments of the distributions is obtained. The corresponding entropy 
indices are calculated indicating chaotic nature of multiparticle production with a 
specific self-similar structure. 

Study of geometrical structure of dynamical fluctuations produced in high^energy col- 
lisions gives an exceptional opportunity to investigate hadroproduction processa The inter- 
mittency/fractal patterns of such density fluctuations observed in all types of collisions, one 
connects with two possible scenarios. A geometrical monofractal structure is associated 13 
with a (most expected) thermal phase transition e.g., from a quark-gluon plasma to nuclear 
matter, while multifractality is assigned^ to a self-similar cascading from partonic state to 
final hadrons with a possible "non-thermal" (non-equilibrium) phase transition during the 
cascade. On the other hand, the quantities under the study are usually calculated through 
averaged moments of the distributions and changes of the density fluctuations from event to 
event are not taken into account. These lead to the loss of information about more structure, 
namely, about chaoticity nature of multiparticle production. Recently, erraticity approach 
has been proposed to study chaotisity in high-energy physicsp in particular, to calculate 
entropy indexes^ large value of which can be considered as a signature of chaotic dynamics 
in self-similar multiparticle production (e.g. in branching processes in QCDtl). 

In this report we perform an investigation of dynamics of fractality and chaoticity in 
multihadron production in central collisions of relativistic nuclei. Such studies have a specific 
interest due to the above mentioned expectation of quark-gluon plasma formation in high- 
energy nuclear reactions and, on the other hand, in a sense, of a possible hadronic nature 
of intermittency, being also stronger-jpth energy decrease™ Note, the study of fractality 
continues our previous investigationsP'H where different regimes in multiparticle productions 
has been mentioned. 

The contribution presented deals with the data came from interactions of the JINR Syn- 
chrophasotron (Duhria) 4.5 A GeV/c 12 C beam with a copper target inside the 2m Streamer 
Chamber SKM-200C3 A central collision trigger was used: absence of charged particles with 
momenta p > 3 GeV/c in a forward cone of 2.4° was required. 

The scanning and the handling of the film data were carried out on special scanning 
tables of the Lebedev Physical Institute (Moscow)LL3 The average measurement error in the 
momentum (e p /p) was about 12%, and that in the polar angle measurements was (&#) ~ 2°. 
In total, 663 events with charged particles in the pseudorapidity window Ar/ = 0.2 — 3.0 
are considered (77 = — lntan($/2)). The accuracy (e v ) does not exceed 0.1. In addition, 
particles with pt > 1 GeV/c are excluded from the investigation as far as no negative 
charged particles were observed with such a transverse momentum. Under the assumption 
of an equal number of positive and negative pions, this cut was applied to eliminate the 
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contribution of protons. The average multiplicity is 23.8 ± 0.4. 

The density fluctuations are considered in the ry-phase-space. To avoid the problem 
connected with a non-flat shape of the distribution p(rj), we use the "cumulative" variable^ 



p(r]')dr}' 



(1) 



with the uniform spectrum p(j]) within the interval [0,1]. 

Self-similar (fractal) dynamics is revealed by a scaling-law, 



Fg cx M v <> 



< <p q < q - 1 , (q > 2) , 
of the qth order normalized ( "vertical" ) factorial momentsl^lO 



F q = 
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over M number of the bins into which the phase space of produced particles is divided. 
Here, n m is the number of particles in the m-th bin per event, N m is the similar number 
but calculated for all M events, and (• ■ ■} denotes averaging over the events, is the 
gth power factorial-multinomial, n(n — 1) • • • (n — q + 1). Note, we apply the modified 
moments proposed c3 to remove the biased estimator of the normalization (especially in 
small bins) , that along with the use of the "transformed" variable (0) allow to study higher- 
order moments. 
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Fig. 1: InF, (H) vs. InM for q = 3 and 5. 



Fig. 2: \ q (W) vs. q. 



The exponents tp q show the strength of an intermittent structure of the distribution. 
Fractality is characterized by the codimcnsions, d q — f q /(q — 1), indicating monofractal 
objects if d q — const. (q), and multifractal ones if the chierarhy, d q > d p , q > p, is satisfied. 
Then, a corresponding connection with a phase transition can be found as discussed in the 
beginning. 

As a signal of the transition, the existence of a minimum of the function 



\ = {Vq + l)/<7 



(4) 



at a certain "critical" value of q = q c is suggestedtl The minimum of eq. (Q) may also be a 
manifestation of coexistence of many small (liquid-type) fluctuations and a few high-density 
onesB 
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Fig. 1 presents dependence of hiF q on InM, depicted for q = 3 and 5. The different 
increase of the moments with increasing M for the different M-regions seeing on the plots 
continues up to q = 8 (not shown) , confirming our earlier results 013 of the existence of 
distinguished regimes of particle production at various bin-averaging scales. 

As in the preceding analysis, multifractality in hadroproduction has been observed, i.e. 

q > p (not shown), supporting a cascading scenario of particle production. 



> 
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Different increase of the d q with q was found, indicating possible change of the regime of 
particle creation. 

Fig. 2 shows the X q function (Q), confirming that at least two regimes of particle pro- 
duction exist: one with the phase transition at q c between 4 and 5, and another one for 
which no critical behavior is reached. The g c -value and the M-intervals, which exhibit the 
minimum o£A g , 11 < M < 17, 11 < M < 24, are found to be about the same as in our previ- 
ous studiesl3'a as well as in recent similar investigations of heavy-ion collisions at ultra-high 
energiestil The most "critical" M-region, 11 < M < 17, is shown hy linear fit in Fig. 1. 

Taking into account multifractality, the critical q c indicates □ a "non-thermal" phase 
transition rather during the cascade than within one phase. Although the interpretation may 
be a matter ijif debate, it must be noted that the minimum was found earlier also in hadronic 
interactionsu at small px and has been indicated in high-energy nuclear interactionsuJ 
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Fig. 3: lnCp, 3 (0) vs. InM. 
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Fig. 4: lnC 2 , 9 (§) vs. In M. 



Intcrmittency has not fully exhausted the charachteristics of dynamical fluctuations. To 
study change of such fluctuations from event to event and then, to find out the chaoticity 
characteristics of the observed self-similar process it was proposeda to consider the normal- 
ized moments C p q , defined by 

C P , q = {f%)/(F q H r , (5) 



of the distributions of the normalized "horizontal" factorial moments 



(f q ). Here, 
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and A is a total number of particles. Note, the order p of the Cp i9 -moments can be non- 
integer, negative or zero, but are restricted with zero f q s ("empty-bin effect"). 

The degree of erratic fluctuations of the self-similar dynamics, indicated by (0), is sug- 
gested to have a scaling behavior, 
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The exponents i/v,? are referred to as erraticity indices and show a strength of "spatial" 
fluctuations. 

The power-law (0) of the C p , 9 -moments vs. M is shown in Figs. 3 and 4. The first one 
illustrates the M-dependence of the moments for different p = 0.4,1.2, and 1.6, and fixed 
q = 3, while Fig. 4 presents the case, when p is fixed, p = 2, and q changes, q — 2,3, 5. As 
seen, the scaling behavior is satisfied up to small bins, however for q > 2 contributions of 
zero /'s (||) are sensitive making the moments to saturate. This also is a reason why we do 
not consider higher g-orders (q > 6) as well as p < 0. 

Note, shown by linear fit the most "critical" M-interval in Figs. 3 and 4 exhibits the 
same rise character as in Fig. 1. Corresponding erraticity indexes are shown in Fig. 5a, 
compared to the "non-critical" fit for 2 < M < 17. The fourth order polynomials are used 
to approximate the ^ Pj9 -slopes as functions of p. The general properties^ of the ^-functions 
are seen: Vp,g < for < p < 1, Vp,<? > for p > 1 (see also Fig. 3). Nevertheless, it 
is clear, that the fits deviates from zero for p = that is connected with the empty-bin 
contributions. 
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Fig. 5: Erraticity indexes Vp-g ([])> 9 = 3,4, as functions of p for different M-intervals: 
(a) 11 < M < 17; (b) 2 < M < 17. 



Of particular interest is an index fi q defined as 

N = d^'«lp=i - (8) 

showniH] to be related to the entropy in event space and be a possible measure of chaoticity. 
The larger n q is, the.more chaotic the system is. In the case of a phase transition high values 
of /x g are expectedl!3 

Table 1 presents the [i q calculated for different M-intervals and q — 2. ..5. The large val- 
ues found for each interval indicate very chaotic dynamics of particle production, confirming 
its cascading nature. It is worthwhile to note increase of the entropy index with approaching 
to the "critical" M-intervals that lends support for a phase transition. However, we must 
emphasize contribution of empty bins at high g's. 
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Table 1. The entropy indexes [i q (g). 



M-interval 


q = 2 


q = 3 


q = 4 


q = 5 


11-17 


0.033 ±0.014 


0.251 ±0.081 


0.926 ±0.218 


2.070 ±0.355 


11-24 


0.032 ±0.012 


0.222 ±0.023 


0.743 ±0.082 


1.581 ±0.170 


2-17 


0.012 ±0.002 


0.095 ±0.010 


0.309 ±0.034 


0.631 ± 0.065 


2-43 


0.024 ±0.001 


0.175 ±0.005 


0.516 ±0.014 


0.967 ±0.025 



In summary, the dynamics of fractality and chaoticity in multiparticle production is 
studied with the pseudorapidity spectra of charged particles produced in central C-Cu colli- 
sions at 4.5 GeV/c/nucleon. The scaled factorial moments are calculated in the transformed 
variables and corrected to take into account the bias of infinite statistics. Multifractal struc- 
ture of pseudorapidity density fluctuations is observed to be of dynamical nature, indicating 
a possible non-thermal phase transition. The existence of two different regimes of parti- 
cle production during the hadronization cascade is mentioned. Xhese results based on the 
increased statistics confirm conclusions of our previous analysisQtl Changes of the dynam- 
ical fluctuations from event to event are further investigated using the distributions of the 
horizontal factorial moments. The normalized moments of these distributions are studied 
and a scaling behavior referred to as erraticity is obtained. The "spatial" entropy indexes 
are calculated shown strong chaotic nature of the hadroproduction process, especially in the 
case of possible non-thermal phase transition. 
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